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^ ' Masses of the ground state light tetraquarks are dynamically calculated in the 

Q ■ framework of the relativistic diquark-antidiquark picture. The internal structure of 

. the diquark is taken into account by calculating the form factor of the diquark-gluon 

CN . interaction in terms of the overlap integral of the diquark wave functions. It is found 
that scalar mesons with masses below 1 GeV: /o(600) (a), Kg (800) (k), /o(980) and 

' ao(980) agree well with the light tetraquark interpretation. 
I ' 

Dh; PACS numbers: 14.40. Cs, 14.40.Ev, 12.39.Ki 

\ The consistent theoretical understanding of the hght meson sector remains an important 

^ I problem already for many years. ^ An extensive analysis of the data on highly excited light 
^ ■ non-strange meson states up to a mass of 2400 MeV collected by Crystal Barrel experiment 
I , at LEAR (CERN) has been published [ij. Classification of these new data requires better 
theoretical description of light meson mass spectra. This is especially important, since 
light exotic states (such as tetraquarks, glueballs, hybrids] within quantum chromodynamics 
(QCD) are expected to have masses in this range 0, 0, Particular interest is focused on 
op \ scalar mesons, their properties and abundance. A generally accepted consistent picture has 
not yet emerged. Experimental and theoretical evidence |3| for the existence of /o(600)((t), 
^ I ii'o(800)(fi;), /o(980) and ao(980) indicates that lightest scalars form a full SU{?>) flavour 
^ ■ nonet. A peculiar feature of their mass spectrum is the inversion of the mass ordering, which 
cannot be naturally understood in the qq picture. This fact stimulated various alternative 
interpretations of light scalars as four quark states (tetraquarks) (1, |6[ in particular diquark- 
antidiquark bound states The proximity of /o/cto to the KK threshold led to the KK 
molecular picture jl]. 

In this paper we calculate the masses of the ground state ((L^)=0) light tetraquarks as 
diquark-antidiquark bound states in the relativistic quark model based on the quasipoten- 
tial approach in quantum chromodynamics. Following 0, 0| the diquark is taken in the 
colour antitriplet state. Recently, in the framework of the same model [ol, 10] we investi- 
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gated the mass spectra of heavy tetraquarks. It was found that many of the newly observed 
charmonium-like states [ll| above open charm threshold, including explicitly exotic charged 
states, could be interpreted as tetraquark states with hidden charm. In the present anal- 
ysis of light tetraquarks we use, as previously, the diquark-antidiquark picture to reduce 



^ A vast literature on the light meson spectroscopy is available. Therefore we mostly refer to the recent 
reviews where the references to earlier review and original papers can be found. 
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a complicated relativistic four-body problem to the subsequent two more simple two-body 
problems. The first step involves the calculation of the masses, wave functions and form 
factors of the diquarks, composed from light quarks. At the second step, a light tetraquark 
is considered to be a bound diquark-antidiquark system. It is important to emphasize that 
we do not consider the diquark as a point particle but explicitly take into account its struc- 
ture given by the form factor of the diquark-gluon interaction in terms of the diquark wave 
functions. 

In the quasipotential approach and diquark-antidiquark picture of tetraquarks the inter- 
action of two quarks in a diquark and the diquark-antidiquark interaction in a tetraquark 
are described by the diquark wave function of the bound quark-quark state and by the 
tetraquark wave function (\I/t) of the bound diquark-antidiquark state, respectively. These 
wave functions satisfy the quasipotential equation of the Schrodinger type [12l] 



where the relativistic reduced mass is 



_ _ - {ml - mlf 



and El, E2 are given by 



= 2M ' = 2M ■ 

Here, M = Ei + E2 is the bound-state mass (diquark or tetraquark), mi^2 are the masses of 
quarks [q = u,d and s) which form the diquark or of the diquark {d) and antidiquark [d') 
which form the light tetraquark (T) , and p is their relative momentum. In the center-of-mass 
system the relative momentum squared on mass shell reads 

, ^ [M' - (mi + m2y][M' - (m^ - m^)'] , 

The kernel l^(p, q; M) in Eq. ([1]) is the quasipotential operator of the quark-quark or 
diquark-antidiquark interaction. It is constructed with the help of the off-mass-shell scatter- 
ing amplitude, projected onto the positive-energy states. For the quark-quark interaction in 
a diquark we use the relation Vgg = Vgg/2 arising under the assumption of an octet structure 
of the interaction from the difference in the qq and qq colour states. An important role in 
this construction is played by the Lorentz structure of the confining interaction. In our anal- 
ysis of mesons, while constructing the quasipotential of the quark-antiquark interaction, we 
assumed that the effective interaction is the sum of the usual one-gluon exchange term and a 
mixture of long-range vector and scalar linear confining potentials, where the vector confin- 
ing potential contains the Pauli term. We use the same conventions for the construction of 
the quark-quark and diquark-antidiquark interactions in the tetraquark. The quasipotential 
is then defined as follows [E]- 

(a) For the quark-quark {qq, qs, ss) interactions, V{p,q; M) reads 



V{p, q; M) = Ui{p)u2{-p)V{p, q; M)ui{q)u2{-q), 



(5) 
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with 



V(p,q;M) = - 



Here, a<j is the QCD couphng constant; D^y is the gluon propagator in the Coulomb gauge, 

47r 



D""(k) 



k2 



(6) 



and k = p — q; 7^ and u{p) are the Dirac matrices and spinors 

u\p) 



e{p) + m 



\ Hp) 



1 

■ p 



X 



(7) 



with e(p) = v^p2 + m^. 

The effective long-range vector vertex of the quark is defined |13] by 



IK 



r^(k) =7^ + — fc = (0,k), 
2m 



(8) 



where k is the Pauli interaction constant characterizing the anomalous chromomagnetic 
moment of quarks. In configuration space the vector and scalar confining potentials in the 
nonrelativistic limit l3] reduce to 



VLAr) = (l-5)Ko„f(r), 



Konf(0 = ^Konf(r) 



with 



Konf (r) = Konf (0 + Konf (0 = At + B , 

where e is the mixing coefficient. 

(b) For the diquark-antidiquark {dd') interaction, V^(p, q; M) is given by 



(9) 
(10) 



nP,q;M) ^ M™2))^«.z.^.(k)i™pil 

+rAP)rAP') U;M^|Krnf(k) + Konf(k)l MQ)^AQ') 



where ((i(P)| J^|(i(Q)) is the vertex of the diquark-gluon interaction which takes into account 
the finite size of the diquark [pC) = ±p) and Q^'^ = (P^o, ±q), Ed = (M^ - M| + 

MJ)/(2M) and E^' = (M^ - Mj + M|)/(2M) " . 

The diquark state in the confining part of the diquark-antidiquark quasipotential ffTTl) is 
described by the wave functions 



i^dip) 



^d{p) 



for a scalar diquark, 

for an axial-vector diquark. 



where the four-vector 

^dij)) = 



[Sd ■ P) 



.£d + 



[^d ■ Pjp 



^d{p)Pt^ = 0, 



(12) 



(13) 
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TABLE I: Masses of light ground state diquarks (in MeV). S and A denotes scalar and axial vector 
diquarks antisymmetric [. . .] and symmetric {. . .} in flavour, respectively. 



Quark 


Diquark 






iviass 






content 


type 


M 


[15] 


[16] 


[17] 


[18] 






our 


NJL 


BSE 


BSE 


Lattice 


[u,d] 


S 


710 


705 


737 


820 


694(22) 


{u,d} 


A 


909 


875 


949 


1020 


806(50) 


[u,s] 


S 


948 


895 


882 


1100 




{u,s} 


A 


1069 


1050 


1050 


1300 




{s,s} 


A 


1203 


1215 


1130 


1440 





is the polarization vector of the axial-vector diquark with momentum p, Ed{p) = y + Mj, 
and £rf(0) = (0, £(i) is the polarization vector in the diquark rest frame. The effective long- 
range vector vertex of the diquark can be presented in the form 

for a scalar diquark, 

(14) 

for an axial-vector diquark. 

Here, the antisymmetric tensor Tj^^ is defined by 

{J:,J^ = ~^{g,,5:-g,^5;), (15) 

and the axial-vector diquark spin Sd is given by {Sd-k)ii = —i^ku] l^d is the total chromo- 
magnetic moment of the axial-vector diquark. 

The constituent quark masses = = 0.33 GeV, = 0.5 GeV and the parameters 
of the linear potential A = 0.18 GeV^ and B = —0.3 GeV, fixed previously p^], have values 
typical in quark models. The value of the mixing coefficient of vector and scalar confining 
potentials e = —1 has been determined from the consideration of charmonium radiative 
decays [l3| and the heavy-quark expansion. The universal Pauli interaction constant k = —1 
has been fixed from the analysis of the fine splitting of heavy quarkonia ^Pj - states |il2|] . In 
this case, the long-range chromomagnetic interaction of quarks vanishes in accordance with 
the flux-tube model. 

At the first step, we take the masses and form factors of the light diquarks from the 
previous consideration of light diquarks in heavy baryons jl^]. The form factor F{r) entering 
the vertex of the diquark-gluon interaction was expressed through the overlap integral of the 
diquark wave functions. Our estimates showed that this form factor can be approximated 
with high accuracy by the expression 

F(r) = l-e-«'"-^'^'. (16) 

The values of the masses and parameters ^ and ( for light scalar diquark [■ ■ ■] and axial 
vector diquark {■ ■ ■} ground states are given in Table [H [TTl 



J, 



2y EdEd 



+ 



2JEdEd 



2Md ^' " 
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TABLE II: Parameters and ^ for ground state light diquarks. 



Quark 


Diquark 




c 


content 


type 


(GeV) 


(GeV2) 


[u,d] 


S 


1.09 


0.185 


{u,d} 


A 


1.185 


0.365 


[u,s] 


S 


1.23 


0.225 


{u,s} 


A 


1.15 


0.325 


{s,s} 


A 


1.13 


0.280 



At the second step, we calculate the masses of light tetraquarks considered as the bound 
states of a light diquark and antidiquark. For the potential of the diquark-antidiquark 
interaction (fTT!) we get in configuration space [lo| 



V{r) 



where 



Vcoui(r) + Vcon{{r) + -< 



+ 



Ei{Ei + M,) E2{E2 + M., 



Mi{Ei + M,] 



1 



M2{E2 + Ma) 
1 



r 



1 1 
2 iMf^Mf 



/V I 
conf \ 



1 



Ei{Ei + Mi) ^2(^2 + Ma 



^Coul(r) 



'L- (S1 + S2) 
1 



/id 








2 






r J 



+- 



EIE2 

1 



Mi(Ei + Mi) M2(E2 + M2) 
L ■ (Si - S2) 



2r 



^Coul(^) + 



4 VMi M2 



L-(Si + S2) 



+- 



+- 



X 



+- 



fid f El E2\ \4o^f\ 



-L -(81-82) 



4 VMi M2 



r + 4 M1M2 Ir ' 



:(8i ■ r)(82 ■ r) - 81 ■ 82 



AV^Coui(r) + 



fil E1E2 



4 M1M2 



AKL(r) 



81 ■ 82 



(17) 



4 Fi(r)F2(r) 

-—as 

3 r 



is the Coulomb-like one-gluon exchange potential which takes into account the finite sizes 
of the diquark and antidiquark through corresponding form factors Fi 2(r). Here, 81^2 and 
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L are the spin operators of diquark and antidiquark and the operator of the relative orbital 
angular momentum. In the following we choose the total chromomagnetic moment of the 
axial- vector diquark /i^ = 0. Such a choice appears to be natural, since the long-range 
chromomagnetic interaction of diquarks proportional to /x^ then also vanishes in accordance 
with the fiux-tube model. 

We substitute the diquark-antidiquark interaction potential (ITTll in the wave equation 
([T]) and solve it numerically in configuration space. The calculated masses of the ground 
states of light tetraquarks considered as light diquark-antidiquark bound systems are given in 
Tables UTT] and HVl In Table UTT] we present masses for light unflavoured tetraquarks (without 
or with hidden strangeness) and in Table IIVI - masses of strange tetraquarks. Possible 
experimental candidates for such states are also given. 

In the diquark-antidiquark picture of tetraquarks both scalar S (antisymmetric in flavour 
[. . .]) and axial vector A (symmetric in flavour {. . .}) diquarks are considered. Therefore we 
get the following structure of the light tetraquark ground (IS*) states (C is defined only for 
neutral self-conjugated states): 

• Two states with J^'" = 0^^: 

X(0++) = SS 
X(0++') = AA 

• Three states with J^'-^ = l"*"^: 

X(l++) 

x(i+-) 
x(i+-') 

• One state with J^'^' = 2++: 

X(2++) = AA. 

The lightest SS scalar (0"*"+) tetraquark states form the SU(3) flavour nonet: one tetraquark 
( [wci] [mc?] ) with neither open or hidden strangeness (electric charge Q = and isospin / = 0); 
four tetraquarks ([sg]['U(i], q = u,d) with open strangeness {Q = 0, ±1, I = \) and 

four tetraquarks ([sg][sg']) with hidden strangeness {Q = 0, ±1, / = 0, 1). Since we neglect 
in our model the mass difference of u and d quarks and electromagnetic interactions, the 
corresponding tetraquarks will be degenerate in mass. 

From Tables IIIII and IIVI we see that the diquark-antidiquark picture can provide a nat- 
ural explanation for the inversion of masses of light scalar 0"^ mesons. Indeed all lightest 
experimentally observed scalar mesons /o(600) (cr), -^"9(800) (/t), /o(980) and ao(980) can 
be interpreted in our model as light tetraquarks composed from a scalar diquark and antidi- 
quark (SS). Therefore, the /o(980) and ao(980) tetraquarks contain, in comparison to the 
qq picture, an additional pair of strange quarks which gives a natural explanation why their 
masses are heavier than the strange Kq{800) [k). Note that physical neutral scalar states 
with 1 = are in fact mixtures of pure tetraquark /o and uo states 

I/) = cos9?|/o) + sin9?|ao). 



V2 



{SA + SA) 



-^{SA - SA) 
-- AA 
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TABLE III: Masses of light unflavored diquark-antidiquark ground state ((L^)=0) tetraquarks (in 
MeV) and possible experimental candidates. S and A denote scalar and axial vector diquarks. 



State 


Diquark 


Theory 




Experiment [19] 






content 


mass 


1 = 


in '^'^ 


1=1 


mass 
















0++ 


SS 


596 /o(600) (a) 


ztnn 1 9nn 




_ 


l+± 


{SA±SA)/y/2 


672 










0++ 


AA 


1 1 7Q 




1 9nn 1 '^00 

L^yjyj -LOyjyj 






1+- 


AA 


1 77*? 
L I 1 O 












AA 


1 

1915 <^ 














J. j.^ f 






{qs){qs) 














0++ 


SS 


992 


/o(980) 


Qson o\ 


nn(QR(W 


984.7(12) 


1++ 


{SA + SA)/V2 


1201 


/i(1285) 






1230(40) 


1+- 


{SA-SA)/V2 


1201 


/ii(1170) 




0\ yLZoo ) 


1229.5(32) 


0++ 


AA 


1480 


/o(1500) 


1505(6) 


ao(1450) 


1474(19) 


1+- 


AA 


1942 


/ii(1965) 


1965(45) 


6i(1960) 


1960(35) 


2++ 


AA 


2097 1 


f /2(2010) 


2011(70) 


f 02(1990) 


2050(45) 


[ /2(2140) 


2141(12) 


\ 02(2080) 


2100(20) 


(ss)(ss) 














0++ 


AA 


2203 


/o(2200) 


2189(13) 






1+- 


AA 


2267 


/ii(2215) 


2215(40) 






2++ 


AA 


2357 


/2(2340) 


2339(60) 







TABLE IV: Masses of strange diquark-antidiquark ground state ((L^)=0) tetraquarks (in MeV) 
and possible experimental candidates. S and A denote scalar and axial vector diquarks. 



State 


Diquark 


Theory 


Experiment [19] 




content 


mass 




mass 


{qq){sq) or {sq){qq) 










0+ 


SS 


730 




672(40) 


1+ 


{SA±SA)/V2 


1057 






0+ 


AA 


1332 


ii:^(1430) 


1425(50) 


1+ 


AA 


1855 






2+ 


AA 


2001 


K2*(1980) 


1973(26) 



\a) = -smv?|/o) + cosv^Ictq), (18) 

where ^ 

fo = —z={[su][su] + [sd][sd]); ctq = [mc?] [mc/] . 
V 2 

However the f — a mixing is small because the Zweig rule is expected to hold in the physical 
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mass spectrum [7|]. Our results support such conclusion, since we get good agreement with 
experimental data already without such a mixing. 

The other scalar tetraquark states can be composed from an axial vector diquark and 
antidiquark (AA). Their masses are predicted to be approximately 600 MeV heavier than 
the SS tetraquarks. The diquark-antidiquark composition also naturally explains the ex- 
perimentally observed proximity of masses of the unflavored ao(1450), /o(1500) and strange 
ii'o(1430) scalars. Let us note that quark-antiquark scalar states are predicted in our model 
to have masses around 1200 MeV {qq) and 1400 MeV {qs). 

The axial vector 1+ states can be composed both from scalar and axial vector diquark 
and antidiquark ((5*^4 ± SA)/^/2) and from an axial vector diquark and antidiquark (^4^4), 
respectively. Our model predicts rather low mass values of the former states composed from 
light quarks ({wc/jf-uJ] ± {ud}[ud]) / \/2, 672 MeV, and of their strange partner ([gsjlwii} ± 
1057 MeV. Such axial vector states are not observed experimentally. Note that 
the recent study [23] also indicates that corresponding light tetraquarks should have masses 
below 1 GeV. On the other hand, there are several candidates for the axial vector ({5's}[gs] ± 
{qs}[us])/V2 tetraquarks both in isospin 7=1 (ai(1260), 6i(1235)) and J = (/i(1285), 
/ii(1170)) channels. However, ordinary qq axial vector mesons are expected to have close 
masses. Therefore the observed states can in principle be mixtures of qq and tetraquark 
states. There are also possible experimental candidates for the axial vector 1~*"~ {qs}{qs} 
tetraquark with isospin 7=1, 6i(1960), and with / = 0, /ii(1965), as well as for the 1^ 
tetraquark, /;,i(2215). 

The diquark-antidiquark ground state tetraquark, composed from an axial vector diquark 
and antidiquark (AA), can also be in the tensor 2+ state. The possible experimental can- 
didates are the following: a {qq}{qq} tetraquark with 1 = 0, /2(1910) or /2(1950); the 
{qs}{qs} tetraquark with 1 = 1, a2(1990) or a2(2080), and with 1 = 0, /2(2010) or /2(2140); 
for the tetraquark /2(2340); for the {qq}{qs} tetraquark i^'2(1980). 

There remains the important problem of describing simultaneously the mass spectrum and 
decay rates of the light and heavy nonets of scalar mesons within the relativistic quark model. 
This requires the inclusion of instanton-induced mixing terms [tI and will be investigated in 
future. 

In summary, we calculated the masses of the ground state light tetraquarks in the diquark- 
antidiquark picture. In distinction with previous phenomenological treatments, we used the 
dynamical approach based on the relativistic quark model. Both diquark and tetraquark 
masses were obtained by numerical solution of the quasipotential wave equations. The 
diquark structure was taken into account by using diquark-gluon form factors in terms of 
diquark wave functions. It is important to emphasize that, in our analysis, we did not 
introduce any free adjustable parameters but used their values fixed from our previous 
considerations of hadron properties. It was found that the lightest scalar mesons /o(600) 
(cr), ^^'0(800) (k), /o(980) and ao(980) can be naturally described in our model as diquark- 
antidiquark bound systems. 
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